import numpy as np
import scipy.linalg as lin


amatrix = [-2.13, 0.51, 0.69,-0.04,-1.05,
            0.81,-1.41,-0.93,-1.40,-0.25,
            1.05,-0.59,-0.66,-1.50,-0.04,
            -1.02, 0.39, 1.32,-1.43, 1.40,
            0.69,-0.52,-0.04,-0.52, 0.68];
A = np.array(amatrix).reshape((5,5)).transpose()


T,Z = lin.schur(A, 'complex')


